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ABSTRACT 

With the imagery from GALEX, HST, 2MASS, and Spitzer, and at the resolution of 
MIPS 24 /xm 6"), we study the variations of the broadband spectral energy dis- 
tributions (SEDs) of star-forming regions within the nearest prototypal major merger 

— the Antennae galaxies. By including MIPS 24 /im dust emission into stellar pop- 
ulation analysis, we reliably, albeit roughly, constrain the star formation histories of 
these 24 /xm selected star-forming regions across the merging disks of the Antennae. 
Our population analysis is consistent with the star formation scenario that, most 
regions across the whole system are at a modest level of star formation with the ex- 
ception of some localized intense starburst sites in the well-known overlap regions 
and the western-loop regions of northern galaxy NGC 4038. Compared with all the 
other regions, the young overlap regions currently (< 10 Myr) are experiencing much 
more violent enhancement of star formation. Across the overlap regions, we suggest 
two sequential star formation paths which we interpret as the imprints of the inter- 
penetrating process of the two merging disks following their second close encounter. 
And we suggest that the star formation in the southern and (especially) northwestern 
edges of the overlap zone may have been just triggered by pre-starburst shocks. The 
well-known mid-infrared ''hotspof in the overlap regions is also a "hotspof at 4.5 /xm, 
whose total 4.5 /xm emission 80% from both hot dust and atomic/molecular lines) 
is comparable with that of the two galactic nuclei. 

Key words: galaxies: individual (the "Antennae" galaxies) - galaxies: interactions 

- galaxies: stellar content - galaxies: starburst - galaxies: photometry 



1 INTRODUCTION 

Galaxy mergers, expecially major mergers, can dramati- 
cally influence the morphological and star-forming prop- 
erties of galaxies over relatively short timescales. Almost 
all the ultraluminous infrared galaxies (ULIRGs) — the 
strongest starbursts i n the local Universe, are in interact- 
ing/merging systems ((Sanders fc Mirabel Il996l ). Moreover, 
IConselice. Chapman &: Windhorst ( 2003 ) suggested that 
about two thirds of submillimeter galaxies at z > 1 are un- 
dergoing major mergers. Galaxy inter actions/mergers seem 
to be very fr equent i n the past (e.g. 'Le Fevre et al.l I2 OO0I: 



Peeters et all 2002: Conselice et al. 2003; Elbaz & Ccsarsk^ 



2003l;'Kartaltcpc et al. 2007^: "de" Rave l et al.ii200a : ILin et al 



'2OO81 : iConselice. Yang fc Bluck i2009i ). Therefore, it is of 
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great importance to understand how the burst of star for- 
mation is triggered in the course of interacting/merging. 

At a distance of 19.2 Mpc (Ho = 75 km s"^ Mpc"^)EI 
, the Antennae (NGC 4038/39, Arp 244) is the nearest pro- 
totypal major merger between two gas-rich s piral galaxies 
(|Toomre fc ToomrjUoT^ i lffibbard et al.ll200ll ). Thus it pro- 
vides us with a unique opportunity to study the induced star 
formation process as a consquence of interaction in detail. It 
has been extensively studi ed at essentially all wavel e ngths 
from X-ray to radio ( iHummel fc van der Hulsg Il986l : 



^ We note the recent debate about the distance to the Antennae. 
ISaviane et all ||2008D determined a distance of ~ 13.3 Mpc from 
the tip of red giant branch, whereas ISchweizer et all (|2008|) esti- 
mated a distance of ~ 22.3 Mpc based on the type la supernovae 
2007sr light curve. Throughout this work, we assumed the tradi- 
tionally adopted Hubble Flow distance. However, our conclusions 
are not affected by the controversy over the distance. 
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(MIR) observations ([Mirabel et al.l 
most intense starburst in this system takes place in the 
so-called overlap region between the two nuclei, which 
indicates its intermediate merging stage, almost totally ob- 
scured in the optical. Considering the abundant molecular 
gas, wide spread star formation and overall modest star 
formation efficiency, Gao et al. (2001) argued that Arp 244 
has the potential of producing an ultraluminous extreme 
starburst in a later stage of merging. 

In the Antennae, observations with HS T have identified 
thou sands of super star clusters (SSCs) (IWhitmore et al] 
I1999I ) possibly being formed as part of the merging pro cess. 
Both theoretical predictions (Goodwin fc Bastian '2006") and 
observations (Whit morc 2004; Fall. Chandar fc Whitmora 
I2OO5I : iMengel eTaLifeoOSh suggest that many of these SSCs 
will dissolve rapidly into the field stellar population in the 
course of galaxy-galaxy merging. This process was coined as 
star cluster's "infant mortalitj/' (|Whitmorell2004l V Further- 
more, a large number of stars may not form in star cluster 
mode. In fact, star clusters found in the deep HST images 
only contribute 9%, 8%, 5%, and 7% of t he apparent total U, 
B, V , and / band light of the Antennae (jWhitmore fc Zhand 
l200d V Due to the uncertainties on the cluster formation his- 
tory and efficiency, recently, iBastian et al.l l|2009l ) pointed 
out the difficulty in the accurate understanding of the age 
distribution of star clusters in mergers like the Antennae. 
Therefore, in order to obtain a complete picture of star for- 
mation histories in the course of merging, we must study 
the extended stellar populations, in addition to these star 
clusters. 

Evolutionary population synthesis has become a pow- 
erful tool of interpretation of the integrated spectropho- 
tometric observations of galaxies and sub-galactic regions. 
The most common method of model-observation compar- 
ison for stellar population analysis in galaxies is SED- 
fitting, with either least-squares or chi-squared minimiza - 
tion technique (e.g. iKong et al.l I2OO0I : ICavazzi et al.ll2002h . 
However, the well-known age-extinction degeneracy prob- 
lem prevents us from obtaining reliable information about 
the star formation history for these galaxies or their sub- 
galactic regions, especially when only broadband photom- 
etry data are available. This is because these extended re- 
gions have an unknown mixture of various stellar popula- 
tions, and the differe nt populations may e xperienc e totally 
different extinctions ( Calzetti. Kinnev fc S torclii-B ergmannI 
ll994l : ICharlot fc FaAoOOl ). 

With the inclusion of the high-quality and high reso- 
lution Spitzer 24 /im dust emission data in our population 
analysis, we show here that the degeneracy between stellar 
population and dust extinction can be broken to a great ex- 
tent. Thus, for the first time, we can reliably, albeit roughly, 
constrain the star formation histories within the Antennae 
galaxies using SEDs over the whole spectral range from far- 
UV (FUV) to MIR. The outline of this paper is the follow- 



ing: In Sect. [2] we introduce the multi-band data that we 
use in this study, and give the multiwavelength photometry 
of star-forming regions selected mainly from the 24 fim im- 
age. Sect. Ogives some brief comparisons of the broadband 
SEDs and their variations across the whole system. Sect. |4] 
presents our methodology to constrain the star formation 
histories across the merging disks, and the main results of 
our population analysis. We discuss these results in Sect. [5] 
and then a summary of our main findings follows in Sect.|Bl 



2 DATA AND PHOTOMETRY 
2.1 Data 

Both FUV (~ 1516 A) and near-UV (NUV; ~ 2267 A) 
images were derived from the G ALEX Ultraviole t Atla s 
of Nearby Galaxies distributed bv ICil de Paz et al.l (|20o3). 
The FWHMs of the PSFs are 5 " and 6" at FUV and NUV, 
respectively. With these data iHibbard et al.l (|2005l ) have 
studied the stellar populations of the famous tidal regions. 

Four broadband (F336W, F439W, F555W, F814W) 
and o ne narrowband (F658N; Ha) images (IWhitmore et al.] 
I1999II taken with the WFPC2 aboard the HST were obtained 
as B associations from the MAST Archive^, and mosaic con- 
taining the four chips for each image was created. 

NIR JHKs atlas images from 2MASS were retrieved 
through the Interactive 2MASS Image Service 0- 

The Spitzer imagery (3.6, 4.5, 5.8, 8.0, 24 /im) of the An- 
tennae was obtained wit h both the Infrared Array Camera 
(IRAC: lFazio et al.ll20oi ) a nd the Muhiband Imaging Pho- 
tometer for Spitzer fMIPS: rRieke et al. 20041 on board the 
Spitzer Space Telescope. The Basic Calibrated Data (BCDs) 
were retrieved with the Leopard software Q . Background 
matching, cosmic-ray removal, flat-fielding and mosaicking 
were performed using the Spitzer Science Center's reduction 
software package MOPEXq Image s of the four I RAC bands 
have previously been presented bv lWang et all l|2004h . 



2.2 Region Selection and Photometry Extraction 

Prior to our multiwavelength photometry comparison and 
extraction, all images were background removed, regis- 
tered/aligned, and resampled to the same pixel scale (1.5"). 
Then all images (except FUV/NUV) are convolved to the 
same resolut ion of MIPS 2 4 ^tm w ith the convolution kernels 
provided bv lCordon et al.l (|2008l ). 

We select star-forming regions primarily as 24 /xm emis- 
sion peaks, since 24 ^m has been shown to be a very good 
local tracer of current star formation (e.g. Calzett i et al.l 
12005' ). Practically, we flrst use the IRAF^ D AOFIND task 
in DAOPHOT stellar photometry package (IStetsonI [l987h 



http:/ /archive. stsci .edu/hst/wfpc2/search. html 

http;/ /irsa. ipac.caltech.edu / applications/2MASS/IM/i nter active. html | 
Availabl e at|http: / / ssc.spitzer.caItech.edu/propkit / spot /| 
http: / / ssc.spitzer.caltech.edu/postbcd / download-mopex.html 
IRAF is distributed by the National Optical Astronomy Obser- 



vatories which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 
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Figure 1. Left: Circular apertures of 9" in diameter were selected mainly as 24 /xm peaks and superposed on color-composite images 
generated from MIPS 24 /xm {red), HST Hq {green), and GALEX FUV(biMe) maps of the Antennae. The HST Ho map has been 
degraded to match the resolution of MIPS 24 ^m. The cross signs mark the nuclear centers of the two colliding galaxies. Pixel scales 
are 1.5". Right: Color composite image generated from the IRAC 8 /^m {red), HST F814W {green) and Ha (blue) maps at their original 
resolutions (re-sampled to 1.5" pixel scales). The pixel scales of the two figures are 1.5". The original resolution of IRAC 8 /^m is ~ 2". 
The mosaicked HST images originally have ~ 0.1" pixels (undersampled). At the adopted distance, 1" corresponds to ~ 93 pc. 



to find all local emission enhancements with 10 a signal-to- 
noise (S/N) ratio threshold in the 24 /im image. Then we ex- 
amine the findings of local enhancements visually and check 
carefully to ensure that these previously detected enhance- 
ments are true star-forming regions rather than artifacts 
due to the Airy diffraction rings or spikes. We also include 
three additional FUV bright regions without any nearby 
24 pm peaks associated (regions 11,12 and 32). Finally, a 
total of 34 nearly non-overlapping circular apertures of 9" 
(~ 800 pc) in diameter (Figure [T| were selected. The large 
aperture size is mainly dictated by the PSF of 24 ^m and 
the rather extended emission features for most regions. 

Our aperture size for photometry is, in any sense, sig- 
nificantly larger than a typical H n region l|Knapenlll998l : 
lOev et al]|2003l ). and even larger than a typical central star- 
forming region of a starburs t galaxy. However, as is already 
found bv lZhang et al.l l|200ll ). the spatial distribution of the 
young star clust ers tends to be corr elated up to physical 
scale of ~ Kpc. iBastian et al. I l|2006l ) have studied several 
star cluster complexes with sizes up to several hundred par- 
sees in the Antennae, they found that the young cluster 
complexes often share the same general velocity distribution 
with associated giant molecular clouds (GMCs), and even 
some complexes themselves are clustered. Our star-forming 
regions are usually spatially resolved into one or few such 
bright star cluster complexes dominated by few bright star 
clusters on the high-resolution optical/N IR imaging (e.g. 
IWhitmore et al.lll999l : iMengel et al.ll2005h . Therefore, most 
of our regions should be supposed to be complexes of star 
clusters over extended background stellar populations. We 



note that some regions show slight, yet not systematic, dis- 
placement (^ 1 — 2") between the 24 /xm and the associated 
FUV peaks. The displacement may indicate that IR and UV 
emission are dominated by different clusters within these re- 
gions. 

Ahhough the measured FWHMs of the GALEX PSFs 
are close to that of the PSFs in MIPS 24 /im image, we 
apply point source aperture corrections to FUV (1.206), 
NUV (1.269) and other bands (1.982). Aperture corrections 
are determined from photometry with increasing aperture 
radii of isolated field stars in the image fields. All the fluxes 
have been corrected for foreground Galactic extinct ion us- 
ing a valu e of E(B- V) = 0.04 6 ( Schlcgcl, Finkbcincr fc Davis! 
1 19981 ) and lCardeUi. Clavton fc Mathij (|1989| ') extinction law 
with R\A=3.1. The uncertainties assigned to the photometric 
values are a quadratic sum of three contributions: variance 
of the local background, photometric calibration uncertainty 
(5% for FUV and WFPC2, 3% for NUV, 10% for IRAC, and 
4% for 24 /im), and poisson statistics noise. 



3 BROADBAND SEDS 

Figure [2] shows the broadband SEDs for some representa- 
tive regions across the Antennae galaxies as compared to 
that of entire Arp 2 20 and the local H n galaxy NGC 2798 
(|Kinnev et al.|[l993l ) . There are some very remarkable differ- 
ences of the SEDs between different regions across the whole 
system of the Antennae galaxies. 

In the IR dust emission bands, 24 /im very small 
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grains (VSGs) and 8/im PAH emission (for high-metallicity 
regions) all show very good correlation with extinction- 
corrected hydrogen recombination emission from starburst 
regions of norm al galaxies to luminou s IR g al axies (LIRGs) 
and ULIRGs (|a1 onso-Herrero et al] l2006l : ICalzetti et alj 
I2OO7I ). Nevertheless, the 8 /xm PAH emission is progres- 
sively depressed relative to 24 /xm with increa s ing star for- 



mation intensities (e.g. ICalzetti et al. 



20051: ISmith et all 
Hence, 



2003) 



the 



I2OO7I : iDraine fc Lil l2007l : Ixhilker et al 
24pm/8/xm ratio would give a rough manifestation of dif- 
ferent strengths of current star formation activities across 
the whole system. 

Obviously, the overlap regions (3 - 9) generally have 
higher 24/^m/8^m ratios than other regions in the An- 
tennae, meaning that the overlap regions are the most 
intense c urrent star-forming sites. Spitzer IRS spectrum 
l|Brandl e t al. 2009) also show that the overlap regions 
have overall the weakest relative strength (normalized to 
15 continuum flux) of PAH features among all these 
star-formin g regions. High-resolution Optical and NI R ob- 
servations jwhitmore et all 1 19991 : iMengel et al.ll2005h sug- 
gest that the overlap regions host most of the youngest clus- 
ters. Region 4, whic h hosts the well-known MIR "hotspof 
l|Mirabel et al.|[l998l ) of the overlap regions, has the highest 
24/im/8/xm ratios, even com parable with that o f Arp 220. 
High-quality MIR spectrum (|Brandl et al.l |2009| ) obtained 
from Spitzer IRS more clearly show that region 4 is 
characterized by very hot dust emission and is among 
the regions with the strongest strength of the radiation 
field. The intrinsically brig htest star cluster (WS95-80, 
IWhitmore fc Schweizer|[l995l ) coincides with region 4. 

The regions (16 - 22) in the northern galaxy NGC 4038 
have the lowest 24/im/8/im ratios and reddest UB colors 
as compared with other regions, which is consistent with 
the fact that large number of rathe r old (100 Myr ^ 500 
Myr) star clusters are found there (|Whitmore et al ] 1 19991 : 
iMengel et ahlbOOSl ) 

In the UV wavebands, on the contrary, the Western- 
loop regions (22-32) have the strongest UV emissions, 
largest UV excesses, and relatively strong 8/im (PAH) emis- 
sions, yet low 24 /im/S fim ratios. These are in contrast 
with those in the overlap regions and indicate their rela- 
tively late star formation stage, wh ich is also evidenced b y 
the large Ha bubbles found there (IWhitmore et al]|l999^ . 
iBrandl et al.l ([20091) also presents the IRS spectra of the UV 
brightest region 26, which has significantly weaker strength 
of the radiation field, yet stronger PAHs emission features, 
than do the overlap regions, in agreement with our broad- 
band analysis. 

The circumnuclear star-forming region of NGC 4038, 
namely region 33, has 24 /Mn ratio close to that of the 

overlap regions, suggesting its intense current star forma- 
tion activity. The prominent, yet redder, UV emission re- 
veals that this region may also hold large post-starburst 
(> 10 Myr) stellar populations. For the two nuclei re- 
gions, NGC 4039 has both weaker 8 ^m PAH emission and 
24 dust emission indicating its lower star formation ac- 
tivity, whereas NGC 4038 has almost the same broadband 
SEDs as that of the local H n galaxy NGC 2798 in the plot- 
ted wavelength range. Thus far there is no strong evidence 
for the presen ce of AGN activit y in the two nuclei. The 
IRS spectrum (jBrandl et al.ll2009l ) suggest that, both nuclei 



regions have significantly weaker strength of the radiation 
field than the overlap regions. In comparison with the south- 
ern nucleus region, the northern nucleus region has flatter 
MIR continuum, stronger PAHs fluxes, yet smaller PAHs 
equivalent widths. These are consistent with our broadband 
analysis, i.e. both nuclei regions have weak current star for- 
mation strength compared to the overlap starburst regions, 
and the southern nucleus region has even weaker current 
star formation strength than does the northern nucleus. The 
stron ger high excitation lines in the southern nucleus re- 
gion (IBrandl et al ] |2009l ) should be attributed to supernovae 
shocks, which ar e also strongly evidenc ed by its very steep 
radio spectrum (|Neff fc Ulvestadll2000l ). We point out that 
our apertures covering the two nuclei (regions 2, 34) are so 
large that the surrounding star-forming regions are also in- 
cluded. 

When we examine the SEDs variations across the over- 
lap regions, we see an interesting phenomenon. From the 
central regions, i.e. regions 6 and 7, to both the north- 
western (e.g. regions 8, 9) and southeastern (e.g. regions 
4) regions, the UV/optical emission gradually become red- 
der and weaker, accompanying with comparable amount of 
energy fraction re-emitted as IR dust emission. In the follow- 
ing sections, we show that the obvious trends of the SEDs 
variations across the overlap regions may reflect the star 
formation modes in the overlap regions. 



4 STELLAR POPULATION ANALYSIS 
4.1 Methodology 

The availability of Ha and MIPS 24 //m, which traces the 
star formation (unobscured and obscured) on timescales 
of ~ 10 Myr, FUV and NUV (tracing star formation 
on timescales of ~ 100 Myr), and the ;7(F336W) and 
5(F439W), whose combination in U/B ratio is sensitive 
(straddles the 4000A break) to the fraction of populations 
younger than a few hundred Myr relative to the old popu- 
lations, makes it reasonable to roughly probe the star for- 
mation histories across the whole merging disks of the An- 
tennae. With minimization technique, we fit the broad- 
band SEDs (FW-Ks) plus the narrowband Ha photome- 
try for the above selected star-forming regions with super- 
postions of three single stellar populations (SSPs), namely 
young (< 10 Myr), intermediate (10 Myr-300 Myr) and old 
(> 300 Myr) populations. In our fitting, we account for the 
different extinctions experienced by different populations. 

For populatio n synthesis models, we use the new ver- 
sion Starburst99 (jVazguez fc Leithere"il 120051 ) SSP models 
with the Padova 2000 stellar evolution tracks, which in- 
clude the full AGB evolution. W e adopt the m ulti-power 
law Kroupa initial mass function (iKroupalbOO^ l. Adoption 
of the bottom-heavy Salpeter IMF would not change our 
main conclusions in this paper . Using the m e tallic ity sen- 
sitive Mg I line at 8806.8 A, IMengel et al] (|2002l ') found 
a handful of star clusters in Arp 244 with solar metallic- 
ity, which is c onsistent with recent metallicity estimation 
(Bastian et al.1 12009) fr om more absorption and emission 
lines in the Antennae. Thus, we flx solar metallicity for 
all the SSPs models. Since we have narrow-band Ha im- 
agery data, and the regions we studied are all actively star- 
forming regions, we include gaseous emission into the SSP 
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Figure 2. Display of FUV-24/jm broadband SEDs of some representative star-forming regions marked in Figure [T] and compared 
to Arp 220 and local H u galaxy NGC 2798. The broadband photometry data points from FUV to iCs, which are fitted with 3-SSP 
superpositions, are represented by filled circles and the Spitzer data points by open circles. For the purpose of co mparison, the broadban d 
WFPC2 data have been transformed to standard C/i?V/ system with the transformation coefficients provided bv lHoltzman et alj l[l99^. 
The best-fit (FVY-Ks) composite stellar broadband SED models have been plotted as solid lines. 




Figure 3. Display of FUV- 24 fira broadband SEDs (extracted from the images with their original resolutions) of four overlap star-forming 
regions, which may have significant starlight contamination by nearby bright regions at short wavebands. Lines and symbols are the 
same as in Figure [2] except that FUV/NUV data are not used in the SED-fittings. It can be seen that the overall shapes of the SEDs 
do not change as compared with the SEDs extracted from the smoothed images (Figure |2j. 



models. Besides the nebular continuum emission and sev- 
eral hydrogen emission lines provided by Starburst99, we 
also account for other strong hydrogen emission lines (i.e. 
Balmer, Paschen, Brackett and Pfund lines) assuming case 
B recombination and strong non-hydrogen element emission 
lines using the line ratios of typical Galact i c H n regions 
compiled by I Anders fc Fritze-v. AlvenslebenI (|2003l ). 

For dust extinction, we use the two-phase du st atten- 
uation recipes developed bv lCharlot fc Fall l|200Cf ). The ef- 
fective absorption curve is proportional to A""'^. One ad- 
justable parameter fi defines the ratio of the total effective 
extinction experienced by intermediate/old (> 10 Myr) to 
young (< 10 Myr) populations. 

It is becoming well known, from the works of the SINGS 
team, that the combination of Ha and 24 /xm accounts for 
both the unobscured and obscured c urrent star formation 
for star-forming regions in galaxies (jCalzetti et al.l l2007l : 



iKennicutt et al.ll2007^ . We adopt the relationship between 
extinction-corrected Ha, observed Ha and 24 fim calibrated 
for the H n regions i n the 33 SINGS galaxies (equation 5 in 



ICalzetti et al. I (120071)') to obtain the Ha extmction, and then 
use the lCharlot fc Falll |200d l extinction curve to derive the 
approximate Ay of the young stellar populations from the 
following equation 



Ay = 2.825 log 



0.031- 



(1) 



where L24^m is the dust-only 24 /xm monochromatic lumi- 
nosity and Lhc is the observed Ha luminosity. 

Therefore, we first restrict the variations of the three 
SSPs to their corresponding age ranges mentioned above 
in order to model the SEDs of various star-forming regions 
across the merging disks in the Antennae galaxies. Then 
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we try to find the best-fit composite models with the 
minimization technique 



X 



E 



ohs.i 



aFrr. 



0'{Fobs,i) 



(2) 



by minimizing Equation [21 In this way, we obtain the four 
best-fit parameters, namely, the mass of the old (Mo), inter- 
mediate (Mi) and young (My) populations, the extinction ra- 
tio (fj,) of intermediate/old to young populations. It should 
be stressed here that we first fix the extinction Ay with 
Equation 1 prior to SED modeling, which afi^ects mostly the 
young stellar populations. This makes our stellar popula- 
tion analysis presented in this paper much less affected by 
the degeneracy between stellar population and extinction. 



4.2 Results of SED-fitting 

4-2.1 Best-fit Parameters 

Our main fitting results are summarized in Table [T] One 
should keep in mind that, we are actually comparing rela- 
tively clustered younger stellar populations with relatively 
extended older populations. Hence, the absolute values for 
the mass ratios between different populations should be sen- 
sitive to the size of photometric apertures and the average 
density of the underlying stellar populations. On that ac- 
count we mainly focus on the relative changes of the mass 
ratios between different populations across the system. 

We notice that the SED-fitting quality for some over- 
lap regions (i.e. regions 8, 9), is not as good as for most 
other regions (Figure (5)). The possible reason can be either 
the significant starlight contamination by nearby bright re- 
gions or the very faint nature (low S/N ratio) of these two 
regions at short wavebands. To check the significance of the 
starlight contamination by nearby bright regions, we also 
present the broadband SEDs for the four optically faint re- 
gions in the overlap region in Figure |3l It can be seen that 
the contamination by nearby bright regions does not change 
the overall SEDs shape for these regions. In fact, when we fit 
the SEDs (excluding the FUV/NUV data) extracted from 
the images at their original resolutions for these regions, our 
main conclusions for the overlap regions do not change. Ta- 
ble [2] lists the fitting results for the overlap regions using 
SEDs extracted from the images at their original resolu- 
tions. According to our best-fitting results for these regions, 
the best-fit ages of intermediate populations < 100 Myr, and 
the best-fitting ages of old stellar population ^ 3 Gyr. 

For the mass ratios of intermediate to old populations 
Mi /Mo , the overlap regions have overall smaller values com- 
pared to other regions, suggesting the overlap star-forming 
regions are very young. Northern edge of the overlap re- 
gions, i.e. region 9, for example, has the lowest Mi/Mo. 
Whereas the western-loop regions of NGC 4038 have overall 
the largest Mi /Mo across the whole system, and they are the 
brighest regions in UV. The western-loop regions host most 
of the intermediate populations. The ratios of most other 
regions across the whole system fall in between the overlap 
and western-loop regions. 

For the mass ratios of young to old populations My /Mo, 
the overlap regions and some western-loop regions have 
My/Mo larger than all the other regions. The similar con- 



trast could also been seen from the age-averaged mass frac- 
tion ratios of young to old populations ( M^'/Agc^ )• over- 
lap and western-loop regions host most of the young popu- 
lations, and they are the most intense current star-forming 
sites across the whole system. 

For the mass ratios of young to intermediate popula- 
tions My /Mi, the overlap regions hold the largest values, 
while the northern regions have overall small My /Mi, con- 
forming with their weak current star formation activities. 
We also point out here one interesting finding. That is, 
across the overlap regions, the northwestern edge, namely 
region 9, and the southeastern edge, region 4, have the 
My/Mi ratios significantly larger than the central regions, 
e.g. region 6. We note that the trend still exist after we 
normalize the mass fraction ratios with their corresponding 
best-fit ages (^^). 

In the last column of Table [1] we also list the 
corresponding total starlight absorbed by dust which 
shoul d be equal to the to tal infrared (TIR) dust emis- 
sion. ICalzetti et al.l ((200^ exploited a relationship be- 
tween 24/im/TIR and 8/24 /im for the star-forming re- 
gions in NGC 5194. Considering the similar global prop- 
erties between the Antennae and NGC 5194, such as si- 
mi liar PA H abundances (the PAH index ~ 4.5 according 
to iDraiiie fc Li (20 07,) du st emission models), comparable 
TIR (ISanders et all 120031) and total molecular gas surface 
densitv l|Wilson et al.|]20o3 ) , we compare the total starlight 
absorption from our models with the TIR estimated from 
dust-onl y 8 ^m and 24 fim using the relation (equation 1 in 
Calzetti et al. 2005) exploited for the star-forming regions of 
NGC 5194. The comparison is shown in Figured Obviously, 
almost all our selected star-forming regions in Arp 244 gen- 
erally follow the same relation derived in the star-forming 
regions of NGC 5194. Considering the general consistency, 
in the following sections, we take the total starlight absorp- 
tion from our SED-fitting models as the best estimate of 
TIR for our selected star-forming regions. 



4.2.2 Star-Forming Regions of Arp 244 on the IRX-UV 
Plane 



I Kong et al.l (|2004l ) found that, while the dust extinction 
is the main factor driving the strong correlation between 
Ltir/Lfuv (hereafter IRX) and UV spectral slope, it is 
the star formation history that afi'ects the degree of devi- 
ation of a star-forming galaxy from the locus of starbursts 
in the IRX-UV color indices' plane. It is surely meaning- 
ful to plot these indices for the star-forming regions of the 
Antennae on the IRX-UV color plane. Figure [S] shows the 
ratio of LTm/ipuv, where TIR is set to be equal to the total 
starlight absorption obtained from our models, as a function 
of /3glx (/3glx = Kong et al. 2004). 

The solid line represents th e tight correlation for starburst 
galaxies (|Kong et al.ll2004l ). Obviously, except the overlap 
regions, almost all the star-forming regions in Arp 244 lie 
below the locus followed by starbursts. The three regions 
that lie above the starburst locus are region 4, 5 and 9. 
Hence, it can be stated here that except for some localized 
starburst sites, mainly in the overlap regions, most regions 
across the whole system are forming stars at a level weaker 
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Table 1. SED-fitting Results. 



Reg. 



fj, Mi/Mo My/Mo My/M; 



My /Age. 
Mj/Age 



■/AgCy 
i/AgOp 



log[Total Abs.](erg/s) 



NGC4039 Arms 



1 


1.1 


0.1 


0.003 


7.0e-4 


0.23 


2.6 


0.5 


41.71 


Southern Nucleus 


2 


2.2 


0.9 


0.03 


6.0e-4 


0.02 


3.8 


1.7 


42.72 


Overlap Regions 


3 


2.5 


0.6 


0.03 


0.003 


0.1 


6.3 


20.2 


43.09 


4 


2.7 


0.3 


0.01 


0.007 


0.7 


10.4 


11.3 


43.21 


5 


4.3 


0.4 


0.002 


0.001 


0.5 


7.5 


12.6 


42.67 


6 


2.3 


0.8 


0.008 


0.002 


0.2 


3.7 


29.0 


42.82 


7 


4.5 


1.0 


0.007 


0.002 


0.3 


8.0 


9.6 


42.53 


8 


4.2 


1.0 


0.001 


7.0(--l 


0.2 


11.1 


6.1 


42.61 


9 


•5.9 


O.iS 


<l.()e-l 


1.0c- 1 


>i.0 >10.0 


:10 


12.40 


Eastern Regions 


10 


3.0 


0.9 


0.01 


5.0e-4 


0.03 


1.5 


5.3 


42.55 


11 


1.8 


1.0 


0.05 


5.0C-4 


0.01 


0.2 


2.1 


42.46 


12 


1.9 


0.9 


0.05 


4.0C-4 


0.008 


0.4 


5.0 


42.55 


13 


2.8 


1.0 


0.01 


7.0e-4 


0.07 


1.3 


3.4 


42.60 


14 


2.1 


0.3 


0.01 


6.0e-4 


0.06 


2.0 


8.7 


42.28 


15 


0.9 


0.5 


0.04 


5.0e-4 


0.01 


1.3 


4.3 


42.36 


Northern Regions 


16 


1.3 


0.1 


0.07 


6.0e-4 


0.008 


1.5 


1.7 


41.57 


17 


1.5 


0.2 


0.1 


6.0C-4 


0.005 


1.1 


1.7 


41.82 


18 


1.1 


0.2 


0.04 


3.0C-4 


0.007 


0.7 


4.3 


41.91 


19 


1.0 


0.3 


0.05 


3.0C-4 


0.006 


0.6 


3.7 


41.89 


20 


1.1 


0.2 


0.07 


4.0C-4 


0.006 


0.6 


3.4 


41.62 


21 


1.2 


0.1 


0.04 


4.0C-4 


0.009 


0.9 


5.8 


41.78 


22 


2.2 


0.1 


0.06 


0.001 


0.02 


1.7 


2.8 


41.92 


WcsLeiii-Looi) Ht'gioiib 


23 


0.5 


0.5 


0.03 


6.0e-4 


0.02 


0.3 


4.4 


42.42 


24 


2.3 


0.7 


0.09 


4.0e-4 


0.004 


0.7 


2.9 


42.09 


25 


3.1 


0.8 


0.03 


8.0e-4 


0.02 


0.8 


11.6 


42.46 


26 


1.9 


0.7 


0.2 


0.007 


0.04 


0.5 


50.0 


42.72 


27 


3.1 


1.0 


0.2 


0.006 


0.03 


0.4 


29.6 


42.72 


28 


3.8 


1.0 


0.2 


0.002 


0.01 


1.3 


17.4 


42.30 


29 


3.1 


1.0 


0.2 


0.003 


0.02 


1.8 


3.0 


42.26 


30 


3.0 


1.0 


0.1 


0.002 


0.02 


0.6 


3.2 


42.43 


31 


2.5 


0.7 


0.03 


7.0e-4 


0.02 


1.2 


7.4 


42.50 


32 


2.6 


0.1 


0.02 


0.002 


0.1 


8.0 


1.5 


42.05 


Circumnuclear Regions of NGC4038 


33 


1.1 


0.8 


0.03 


4.0e-4 


0.01 


0.3 


1.3 


42.90 


Northern Nucleus 


34 


1.4 


0.6 


0.01 


7.0e-4 


0.05 


1.5 


1.4 


43.03 



Xr is the reduced minimum value, fi is the extinction ratio of intermediate/old to young populations. My . 



Mi and Mo represent the mass of young, intermediate and old populations respectively. 



My/Agey 



and 



My/AgOo 



Mi/Agei Mi/Ageo 

represent the corresponding age-averaged mass fraction ratios. The total starlight absorption obtained from 
the models is listed in the last column. 



Table 2. Fitting Results for the SEDs (without FUV/NUV) extracted from the images at their original 
resolutions. Only results for some of the overlap regions which may have significant mutual light contamination 
are listed. 

Overlap Regions 



Reg. x'r M Mi/Mo My/Mo My/M; mJ/ZZ l°g[Total Abs.](erg/s) 



3 


1.0 


0.5 


0.04 


0.01 


0.3 


5.0 


40.0 


42.97 


4 


0.3 


0.3 


0.01 


0.05 


4.5 


9.4 


100.1 


43.50 


5 


0.7 


0.8 


0.004 


0.005 


1.3 


4.7 


13.5 


43.11 


6 


1.9 


1.0 


0.05 


0.01 


0.2 


1.5 


73.0 


42.94 


7 


1.8 


0.3 


0.02 


0.004 


0.2 


2.9 


42.5 


42.54 


8 


0.9 


1.0 


0.04 


0.002 


0.05 


0.6 


8.8 


42.82 


9 


1.0 


0.7 


0.002 


0.001 


0.5 


5.0 


4.2 


42.52 
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Figure 4. Total starlight absorption obtained from the SED- 
fitting models is compared with total infrared luminosity esti- 
mated from dust-only 24 [im, and 8 [im, usi ng the relation cali- 
brated in NGC 5194 bv lCalzetti et al.l Hooi) (hereafter TIRcOs). 
The line with slope of unity and null intercept is also plotted 
{solid line). 



than that of the starbursts. This is also consistent with our 
further analysis on star formation histories. 

Figure [6] plots the ratios of L{Hoicorr) / L{Ks) vs. 
Lx{U)/Lx{B). Star-formin g regions in NGC 51 94 selected 
(mainly as 24 /xm peaks) bv lCalzetti et al] l|2005h are shown 
as small pluses. Also plotted is the archetypical starburst 
galaxy NGC 7714 as a whole {large cross). The modest to- 
tal star formation rate {^-^ 3.4 M0yr~^) and star forma- 
tion intensity (~ 0.015 Mgyr-^Kpc"^) of NGC 5194 place 
it among the quiescently star- forming galaxie s , alth ough it 
hosts a LlNER-type nucleus. ^Calz etti et al. I (|2005l ) found 
that the star-forming regions in NGC 5194 have properties 
quite similar to that of the normal star-forming galaxies, 
rather than that of starbursts. Assuming L\{U) / Lx{B) and 
L{Hacorr) / L{Ks) roughly represent the ratios of recent-to- 
past and current-to-past star formation strengths, respec- 
tively. Figure |6] reveals similar results as that of the IRX- 
UV diagram. Namely, almost all the star forming regions of 
Arp 244 have current star formation strength comparable 
with star-forming regions in quiescently star-forming galax- 
ies (low L{Hacorr)/ L{Ks)), except for some of the overlap 
and western-loop regions. 



4-2.3 Emission Excess at IRAC 3.6^.m and 4.5iJ.m 

Since we have done relatively sophisticated SED-fitting for 
these star- forming regions, we could estimate the hot dust 
emission at IRAC 3.6 /im and 4.5 ^m which are often as- 
sumed to be dominated by stellar photospheric emission in 
literatures. With the average of 3.6 /xm and 4.5 ^.m fl u xes as 
the underlying stellar continuum level, IWang et all (|2004l ) 
concluded that most of the emission (>90%) in these two 
bands comes from direct stellar contribution for the whole 
system. From our stellar SED-fitting, we could easily see 



O 



1 r 



"^Overlop Regions 
"^Western — loop Regions 
■ 4038 Nucleus 
□ 4039 Nucleus 
A East Regions 
ANorth Regions 
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#Circum 4038 




-2 
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Figure 5. Ratio of TIR (set to total starlight absorption from 
models) to FUV luminosity as a function of equivalent ultraviolet 
spectral slope (the IRX-UV diagram). Data for different star- 
forming regions within the Antennae are plotted with different 
symbols. The solid line resprese nts the mean IRX -UV relation 
for 50 starburst galaxies fitted bv lKong et al. I 1I2OO4I') . 



(Figures [2] and [Sjl the prominent emission excess at 3.6 /im ( 
~ 10% - 65%) and 4.5 fim { ~ 10% - 80%) for most of these 
star-forming regions relative to the stellar (plus nebular con- 
tinuum) emission, especially for the overlap regions. The sig- 
nificant emission excess at 3.6 fim and 4.5 /xm brings out the 
caution needed in using them as representation for the un- 
derlying stellar mass, particularly when studying the active 
star-forming regions. Since we have largely accounted for hy- 
drogen/helium recombination emission in our SED-fitting, 
according to the existing observatio ns of spectra in these two 
band s of compact H n reg ions (e.g. lMartin-Hernandez et alj 
I2OOOI ; iPeeters et all 120021 ). the other possible main mecha- 
nisms of the emission excess could be vibrationally excited 
molecular hydrogen emission, atomic fine-structure lines, or 
very small hot dust (e.g. 3.3 fim PAH feature) emission. 

For region 4, the MIR ''hotspof , the excess emission 
components contribute 65% and 80% of the total emis- 
sion at 3.6 ^.m and 4 . 5 ^,ra respectively. As is already pointed 
out by IWang et al.l (|2004l '). the extremely high obscuration 
of region 4 must be responsible for the large excess to 
some extent. The high-quality MIR spectrum of this re- 
gion obtained from IRS on board Spitzer is ch aracterized by 
prominent PAHs and very hot VSGs emission (|Brandl et alj 
I2OO9I ). Given the prominent PAH feature emission of re- 
gion 4, the exceptionally low excess emission fiux ratio of 
3.6 /im/4.5 /im (~ 0.9) can not be accounted f or by the 
pure dust emission models of iDraine fc Further- 
more, the K-hand spectra l|Gilbert et al. 200oh in the obser- 



^ Considering the hot dust emission contribution to Ks band, our 
estimation should be lower limits. 
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Figure 6. Ratio of extinction-corrected (using 24 /ivci ) 
L{Hacorr) to L(A's) vs. log[Lx{U)/L^{B)]. Th e symbols are 
as Fig ure [5] The star-forming regions selected by ICalzetti et al.l 
(|2005|) in the quiescently star-forming galaxy NGC 5194 are 
shown by small pluses. The archetypal starburst galaxy NGC 7714 
as a whole is represented by large cross. Obviously, just like the 
majority of star-forming regions in NGC 5194, most regions (ex- 
cept some of the overlap and western-loop regions) in Arp 244 
host a significant fraction of older stellar populations and have 
relatively lower current star formation strength compared to typ- 
ical starbursts. See the text for further explanation. 



vations of the compact star clusters coincident with the MIR 
"hotspof are characterized by prominent nebular and fluo- 
rescent H2 emission with slightly rising continuum toward 
the red (due to hot dust emission). The fact that region 4 has 
very flat radio spectrum (Neff fc Ulvestad.2000. ) and promi- 
nent PAH feature emission excludes the shock origin of the 
possible molecular/atomic emission. Thus, the fine-structure 
lines from heavy elements, the molecular lines from FUV fiu- 
orescence in dense photodissociation regions (PDRs) and hot 
dust (including PAHs) emission may all have significant con- 
tribution to the emission excess at IRAC 3.6 fim and 4.5 /im. 

Interestingly, the total 4.5 /im flux of region 4 is even 
^ 20% higher than that of the southern nucleus and nearly 
equal to that of northern nucleus. In short, the excess emis- 
sion (either from hot dust or molecule/ionized atoms) at 
IRAC 4.5 /xm bands have nonneglig ible (> 18% ) contri- 
bution to the total flux for the whole system. 



* This is obtained by the sum of flux of excess emission for all 
the star-forming regions selected in this work divided by the total 
fiux for the whole system. So here our estimation is a lower limit 
for the whole system. 



5 DISCUSSION 



5.1 Star Formation Histories Across the Whole System 

Since we have the mass ratios between the three stellar popu- 
lations for these star-forming regions, using the best-fit ages 
of the three populations for each region, we could easily de- 
rive the ratios between the age-averaged star formation rate 
(SFR) for the three (i.e. current, recent and past) star forma- 
tion epoches related to our three populations (Table [1] and 
2). We find that, except for some of the overlap regions and 
the western-loop regions of NGC 4038, which have the ratios 
of current-to-past age-averaged SFR (defined approximately 
10, almost all the other star-forming regions 



as '«~c„i~ 



Mo/AgCo 

have the ratios > 1. Specifically, the regions which have the 
ratios > 10 are 3, 4, 5, 6, 25, 26, 27, and 28. Hence, it can 
be said that, some of the overlap and western-loop regions 
are experiencing bursts of star formation, whereas most of 
the other regions are forming stars at a moderate level com- 
parable with normal star-forming galaxies. The location of 
these star- forming regions on the IRX-UV plane (Figure[5]) is 
also consistent with the fact that most of the regions are at a 
modest star formation intensity. For the ratios of current-to- 
recent age-averaged SFR (defined as ^^yxf^) Table [T|, 
the overlap regions have the ratios > 3, whereas almost all 



My/AgCy 

Mi /Age; 



the other regions have the ratios ~ 1. The high 
for the two nuclei may be attributed to their surrounding 
star-forming regions. Hence, comparing with the overlap re- 
gions, the star formation of most regions across the whole 
system is not only modest (except the western-loop regions) 
but also continuous during the recent tens of million years 
(typical best-fit ages of the intermediate stellar populations 
for these regions). Namely, while both the overlap regions, 
and the western-loop regions are the most intense current 
star-forming sites across the whole system, the overlap re- 
gions are now experiencing much more violent enhancement 
of star formation compared to all the other regions. 

In fact, the star formation histories for these star- 
forming regions could also be simply probed through the ra- 
tios of L\{U) / L\{B) and L{Hacorr)/L{Ks), which roughly 
represent ratios of the recent-to-past star formation strength 
and the current-to-past star formation strength, respec- 
tively. Figure |B]clearly shows that, except for the overlap and 
western-loop regions, most regions are forming stars at mod- 
erate level (low L{Hoicorr) / L{Ks)) comparable with star- 
forming regions in quiescently star-forming galaxies. And 
meanwhile these regions host a remarkable fraction of ag- 
ing stellar populations (small L\{U)/Lx{B)), just like the 
star-forming regions in NGC 5194. The overall larger ra- 
tios L\{U) / Lx{B) for the overlap, the western-loop and the 
eastern (to a lesser extent) regions compared with the star- 
forming regions in quiescently star-forming galaxies indicate 
that these regions may have just experienced a period of in- 
tense starburst in the recent few hundred million years. It is 
also notable of the much higher ratios of current-to-past star 
formati on strength f o r som e of the overlap and western-loop 
regions. IWang et al.l l|2004l ) also got similar results from the 
flux ratio of the dust-only IRAC S/xm and the underlying 
stellar continuum. 
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5.2 The 20 cm-to-CO ratio map as star formation 
efficiency map 

Taking advantaR e of the tight correlation between FIR an d 
radio continuum (|Condonlll99l lYun et al.|[200ll : [Beiill2003D , 

which ap pears to be va hd at least on kilo parsec scales in 
galaxies (|Lu et al.|[l99^ . iGao et al.l l|200ll ) constructed a 
star formation efficiency map using the 20 cm-to-CO ra- 
tio. Overall, our results of population analysis are consistent 
with their star formation efficiency map, viz, except some 
localized starburst sites mainly in the overlap and western- 
loop regios, most regions across the whole system are form- 
ing stars at a quite moderate level comparable with normal 
star-forming galaxies. This verifies the practice of using ra- 
dio continuum as an indicator of the star formation and the 
radio-to-CO ratio maps as representation of star formation 
efficiency maps in most cases. 

Nevertheless, some exceptions do exist, like for some 
overlap regions of galaxy pairs in high-speed collision. For in- 
stance, in the Taffy galaxy a large portion of the synchrotron 
radio emission may be related to gas collision shocks, rather 
than supernovae remnants (S NRs, most likely related to 
recen t star formation) shocks (jCao et al.l l2003l : IZhu et al 
120071) . In the overlap regions of the Antennae, iGao et al 
||2001I 1 found that the radio-to-CO ratios progressively in- 
crease from the southeastern side to the northwestern edge 
across the overlap regions. However, we note that, the lower 
radio-to-CO ratio of the southeastern region (4) compared 
to other overlap regions may be, rather than due to lower 
star formation efficiency, assigned to other various causes. 
First, the very flat radio continuum suggests that the 20 cm 
continuum in region 4 , rather than being dominated by 
synchrotron emission from (maybe) SNRs shocks like other 
overlap regions, is primarily thermal free-free emission from 
young, compact H n regions. This means that, in such vio- 
lent overlap starburst envirionment, large numbers of super- 
novae events associated with current star formation epoch 
have not happened. Second, studies o n the local IR/radio 
correlation within n earby galaxies (e.g. iHughes et al.|[200^ : 
iMurphv et al.|[200^ ') demonstrate the weak trend of increas- 
ing IR/radio ratio with increasing IR luminosity within in- 
dividual galaxy. Both our TIR estimates (equal to the total 
starlight absorption) and recent I R estimate b ased on the 
15 iim and 30 pim continuum fluxes l|Brandl et al.,. 2009. 1 sug- 
gest that region 4 has the largest IR luminosity across the 
whole system. Finally, the very hot dust emission, as we have 
shown above, for the southeastern side (i.e. region 4) implies 
that most energy there may emit in the MIR regime, which 
again is different from most of the other regions. Therefore, 
unlike the usually continuous star formation in a relatively 
long timespan (^ 10* yr), in the violent galaxies interaction 
regions, both the strong shocks as a result of cloud-cloud 
collisions and the strong variations of star formation rate 
in short timescales (tens of Myr) may all make the 20cm- 
to-CO ratio fails to be an efficient star formation efficiency 
indicator. 



5.3 Sequential Star Formation Paths in the Overlap 
Regions 

Across the overlap regions (Table[T]and 2), the northwestern 
edge and the southeastern edge have both higher mass ratios 



of young to intermediate populations My /Mi (i.e. regions 9 
and 4) and higher correspondingly age-averaged SFR ratios 
^'^/Agc'' than do the central regions (e.g. regions 6). We 
like to interpret the trends as two sequential star formation 
paths. One is from the central regions (i.e. regions 6) to 
the southern edge (e.g. region 4), and the other is from the 
central to the northwestern edge (e.g. region 9). 

To check if there are any trends for the spatial distribu- 
tions of the star clusters across the overlap re gions, we re- 
fer to the cluster distribution maps derived bv lZhang et al.l 
(!200.l|). We find a slightly proportionally deficit of clusters 
with ages of tens of Myr for both the northwestern (corre- 
sponding to our region 9) and the southern ed ges compared 
to the central regions (e.g. region 6). Recently. iMengel et al.l 
(j2005l ) also obtained the spatial distributions for clusters 
with Br7-determined and CO-index-determined ages. Inter- 
estingly, the distributions also show the relative deficit of 
intermediate-age (~ 10 Myr) clusters for both the north- 
western and southern edges, although the sample size and 
the probed age ra nge are all small compared to that of 
IZhang et al.l (|200ll ). The cluster distributions are in agree- 
ment with the sequential star formation paths mentioned 
above, i.e. both the northwestern and southern regions are 
just beginning their recent star formation episode. 

|jog fc SolomonI l|l992l ) proposed a physical mechanism 
to explain the origin of the enhanced star formation occur- 
ring in situ in the overlapping regions of a pair of colliding 
galaxies like the Antennae. In their model, following a col- 
lision between galaxies, the H I cloud-cloud collisions from 
the two galaxies lead to the formation of hot, ionized, high- 
pressure remnant gas that compresses the outer layers of pre- 
existing GMCs in the overlapping regions. This makes the 
GMCs sheUs become gravitationally unstable, which trig- 
gers a starburst in the initially barely stable GMCs. Al- 
though generally H I cloud-cloud collisons should be more 
efficient than that of the GMCs due to the much smaller 
mean free path of an H I cloud than CMC, the huge con- 
centration of molecular gas in the ove rlap regions of Arp 244 
j Wilson et al.ll2000l : ICao et al.|[200lh may make the GMCs 
collisions also possible. 

Based upon this model, across the overlap regions of 
Arp 244, the sequential star formation paths we have found 
can be explained naturally. As the two colliding disks begin 
to interpenetrate each other, the regions that overlap first, 
e.g. region 6, may have more of the gravitationally unstable 
GMCs layers formed first, leading to in situ starburst first 
there. As the colliding/merging proceeds and the overlap- 
ping zone between the two colliding disks bulks up, cloud 
collisions from the two colliding disks spread to more and 
more regions, i.e. both the northwest (e.g. to regions 8, 9) 
and southeast (e.g. to regions 5, 4) of the overlap regions. 
This leads to progressively lagging starbursts triggered by 
the radiative shock compressions toward both the northwest 
and southeast directions of the overlap regions. In short, the 
identified sequential star formation paths could be the im- 
prints of the interpenetrating process of the two colliding 

galaxy disks. 

Both kinematic analysis (iHibbard et al.ll200ll) and nu - 



merical simulations (|Toomre fc Toomre 11973 : lBarnej[l988h 



suggest that the two galaxies of the Antennae system be- 
gun their first close encounter several hundred million years 
ago, and the two colliding galaxies may have passed their 
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first close encounter l|Mihos. Bothun fc" Richstond Il993h . 
The first peak of larse-sc ale starburst phase may Ii ave just 
passed in the Antennae ijMihos fc Hernauistlll99a ). which 
can be evidenced by the moderate star formation strengths 
shown in Figs. 5 & 6, consistent with t he star formatio n 
efficiency map across the whole system (jCao et al.ll200lt ). 
While the moderate, continuous (recent) star formation for 
most regions is consistent with this scenario, the currently 
violent enhancement of star formation and the sequential 
star formation paths for the overlap regions may suggest 
that now the two colliding galaxies are just launching their 
second close encounter. 

After analysing H2 f = — S( S) X = 9.66 nm line 
emission obtained by ISOCAM CVF, iHaas. Chini fc Klaad 
found that both the southwestern and northwestern 
edges of the overlap regions, which are very close to the star- 
forming region 5 and regions 9, respectively, have exception- 
ally high L{H2) / L{F I R) ratios that exceed that of all other 
known galaxies. But the absolute current /recent star forma- 
tion there (especially the northwestern edge) are very weak, 
just as we find from population analysis (see Table [T] and 
2). They suggest that the high H2 emission there should be 
excited by pre-starburst shocks which are caused by cloud- 
cloud collisions. The low mass fractions of intermediate pop- 
ulations and the very high (age-averaged) mass ratios of 
young to intermediate populations for the two edges, espe- 
cially the northwestern edge, indicate there indeed are very 
young star-forming sites, which are most probablely to be 
triggered by pre-starburst shocks following the second close 
encounter between the two galaxie s. However, recent h igh- 
quality observations from Spitzer (|Brandl et al.l |2009| ) de- 
tected about five times less integrated H2 S(3) line flux than 
does ISOCAM CVF, and didn't find the previously claimed 
strong H2 emission peak in the northern overlap zone. These 
new observations cast doubt on the pre-starburst shock ori- 
gin of the H2 emission in the overlap regions. 



6 SUMMARY 

To summarise, taking advantage of the availability of 
multiwavelengh imagery from FUV to 24 ^m from 
GALEX, HST, 2MASS and Spitzer in both high resolution 
and high sensitivity: 

• We compare the broadband SEDs of star-forming re- 
gions selected as 24 /im peaks across the whole merging disks 
of the Antennae galaxies, which provides us a basis to com- 
prehend the complete picture of star formation histories. 
The large ratios of 24 /im/8 /im for the overlap regions and 
the blue, strong UV emission for the western-loop regions 
demonstrate that currently they are the most intense star- 
forming sites, although the western-loop regions are at a rel- 
atively later star formation stage compared to the overlap 
regions. Most of the other regions, which have redder UB 
color, weaker IR dust emission and UV emission, across the 
whole system are forming stars at a quite moderate level 
during the past ~ 100 Myr. 

• We roughly constrain the star formation histories of 
these active star-forming regions across the whole system, 
with the degeneracy between stellar population and extinc- 
tion broken, by including 24 /xm dust emission into pop- 
ulation analysis. Compared with other regions, the over- 



lap regions are now experiencing much more violent en- 
hancement of star formation, although both the overlap and 
the western-loop regions are the most intense current star- 
forming sites across the whole system. Our analy si s is i n 
general agreement with the findings of ICao et all l|200l[ ) , 
i.e. except for some localized violent starbursts confined 
mainly in the overlap regions and the western-loop regions of 
NGC 4038, the bulk of star formation is at a moderate level 
comparable to that of star-forming regions in quiescently 
star- forming galaxies. 

• We suggest two sequential star formation paths across 
the famous overlap regions, which may reflect the (second) 
interpenetrating process of the second passage between the 
two colliding galaxy disks. We also suggest that the recent 
star formation of both the northern and southern edges of 
the overlap zone might be just triggered by pre-starburst 
shocks. 

• We report the nonnegligible (> 18%) excess emission 
contribution to the total IRAC 4.5 /xm for the whole sys- 
tem. The well-known brightest MIR "hotspof in the overlap 
regions has total 4.5 /xm emission (^ 80% excess emission) 
even higher (by ~ 20%) than that of the nuclear region 
of NGC 4039 and nearly equals to that of the nuclear re- 
gion of northern galaxy NGC 4038. The unusually low ra- 
tio of 3.6 /im/4.5 fj,m implies that, in addition to hot dust 
emission, other emission mechanisms, such as atomic flne- 
structure lines and vibrationally excited molecular hydro- 
gen lines from dense PDRs, have signiflcant contribution to 
IRAC 4.5 Atm. 
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